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System components and tunable parameters , 



Program 


Function performed 


BlastParse.pl 


Parse BLAST output 


mark substitution__heteroplasmy.pl 


Mark ambiguous base calls 


erfract_SE_consensus,pl 


Extract a consensus sequence from a 
phrap polyphred run. replacing low quality 
bases with and polyphred rank 1 
bases with 'N\ 


compute_coverage.pl 


Calculates overall, "forward strand, and 
reverse strand coverage for all bases in 
the consensus sequence 


count_hv1_deletes.pl 


Count deletions in HV1 relative to the 
rCRS 


border_index.pt 


Compute start and end positions of HV1 
and HV2 regions. 


fix_mitotype_.reporting_range.pt 


Output a mtDNA profile based on input 
variants list and reportable range. 


find _badjraces _from_blast_report.pl 


Reject reads that do not align 
appropriately to the rCRS 


determineReadTypes.pl 


Adds template name, template type, and 
primer type to phred output files. 


seq2delta_vs 


Align mtONA profile to rCRS and report 
variants according to nomenclature. 


calculate_coverage_mitotype.pl 


Mask mtDNA consensus sequence based 
on required coverage and calculated 
coverage 


count_hv1_inserts.pl 


Count insertions in HV1 relative to the 
rCRS 


generate_hv_mask_fasta_files.pl 


Extract the HV1 portion and HV2 portion 
of the mtDNA consensus sequence, 
based on the computed HV1 and HV2 
regions. 


flipjasta.pl 


Invert an X masked sequence output by 
cross match 



FIG. 5 
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'tVedo&iiustcbase coraclly called 

Differences take into account the consensus of the two analysts 



Total Thie Positives: = N^^i 


.46 


10 


Total True Negatives: A WXOfflatrf = A,^, 


77,358 


110,354 


Total False Positives: N 4UIfle ^= A**^ 


95 


807 


Total False Negatives: A^ — i^N, — j • 


6 1 


14 


Total Incorrect 


0 


0 


Sensitivity: TP/(TP+FN) 


7173% 


41.67% 


•Specificity: TN/fTN+FP) 


99.88% 


99.27% ' 


Positive Predictive Value: TP/(TP+FP) 


14.41% 


. 1.22% 


Negative Predictive Value: TN/(TN+FN) 

*• 


99,99% 


99.99% 



FIG, 6 
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